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ABSTRACT 

The paper investigates the overall and detailed features of cosmic ray (CR) 
spectra in the knee region using the scenario of nuclei-photon interactions around 
the acceleration sources. Young supernova remnants can be the physical realities 
of such kind of CR acceleration sites. The results show that the model can well 
explain the following problems simultaneously with one set of source parameters: 
the knee of CR spectra and the sharpness of the knee, the detailed irregular 
structures of CR spectra, the so-called "component B" of Galactic CRs, and the 
electron/positron excesses reported by recent observations. The coherent expla- 
nation serves as evidence that at least a portion of CRs might be accelerated at 
the sources similar to young supernova remnants, and one set of source param- 
eters indicates that this portion mainly comes from standard sources or from a 
single source. 

Subject headings: cosmic rays - knee - fine structures - e"''e~ excesses- "component 
B" -single source 



1. introduction 



Since the discovery by lKulikov fc KristiansenI ( 1l958l ). the "knee" of cosmic ray spectra 
has been one fundamental problem of CR physics for half a century. Many theoretical works 
try to explain this interesting and important phenomenon. The most popular explanation 
attributes the knee to th e inefficient acceleration of the Galactic CRs by the accelerators 
above PeV energies (e.g., iKobayakawa et al.ll2002l : ISveshnikoval l2003t lErlykin fc Wolfendale 



20011). Alternative possibilities inc l ude the leakage of C Rs when propagating in the Galaxy 



f lPtuskin et al.lll993l: iRoulet 



2004 



background light (IKarakula fc TkaczykI Il993l : ICandia et al.l l2002l ) or neutrinos (jWigmans 



Lagutin et al.l 120011). interactio ns between CR s and the 
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20031 ) before arriving at the earth, or exotic interaction of C Rs in the atmosphere where 
undetectable particles are p roduced and missing the detection ( iNikolsky fc RomachinI 12000 
Kazanas fc NicolaidisI 120011 ) . 



From the experimental aspects, the measurements of the CR spectra around the knee 
region become increacingly precise, which can even reveal some fine structures of the knee. 
After a long term operation, the Tibet Air Shower array reported a very go od measurement of 



the k nee spectra with unprecedented high statistics and low systematics (lAmenomori et al. 
20081 ). Especially interesting signature of the Tibet result is that the CR spectra show 



a very sharp break of the spectrum index around 4 PeV. At almost the same time sev- 
eral experir aents have report ed their new measure ments with the simila r behavior , such as 
KASKADE jApel et al.ll2009[). ARAGATS-GAMM A JCarvaka et al.l boosh. Yakutsk Jivanov et al, 



20091 ). and MAKET-ANI flChilingarian et a: 



tional interpretations of the knee (jShibata 



2007 ). Such a sharp knee challenges the tr adi- 
20091 lErlvkin fc Wolfendal] boooi IhuI|2009I ). It 



is shown that if adopting an exponential-like cutoff of each component of CR specie s with 
low He flux, i t w il l be v ery difficult to reproduce the sharp knee data ( 1Shibatall2009l ). It is 
suggested by IHuI (120091 ) that a double power-law may well fit the observational data with 
high He fiux, which indicates that He may be the main component around the knee. Fur- 
thermore Hillas suggested that there should be another Galactic component, "component 
B", to explain the CR spectra above 10 PeV JffiUaslboOsh . 



Besides the sharp transition of the knee. lErlykin fc Wolfendald ( 120091 ) carefully analyzed 
the CR spectra of individual experiment. By renormalizing the energy with respect to the 
break point (measured knee energy)of individual experiment, the problem related to the 
uncertainty of the absolute energy scale can be avoided, so the deviations of the observed 
spectra from the fitted spectra can be combined for all experiments. The result clearly 
shows the peculiarities at the positions expected for CNO group and Fe group if the knee 
is corresponding to the position of He. Interestingly the energies of these fine bumps are 
proportional to the mass number of the several major nuclei species: proton. He, CNO and 
Fe. The sharp knee and th e irregularities of CR spect ra are regarded as evidence for the 
single source origin of CRs (lErlykin fc Wolfendald l2009[ ). 



Another important developm ent in CR physics is the new discovery of electron /positron 



Adriani et al. 


2009: 


Chane et al. 


2008 





20081 . 12OO9I : lAbdo et al.l 120091 ). To explain the positron fraction and electron spect rum ex- 
cesse s simultane o usly o ne may need to introduce some exotic sources of e"'"e~ pairs (ISerpico 
2OO9I ). IHu et al.l (I2OO9I ) (hereafter Paper I) proposed a model resorting to e"'"e~ pair through 
interactions of CR nuclei and ambient photons around the acceleration sources, which can 
explain the knee of the CR spectra and the electron/positron excesses at the same time. 



- 3 - 



Based on that model we further study the detailed structures of the CR spectra in this work, 
intending to reproduce the sharp knee and fine structures mentioned above. In our inter- 
action model the threshold energy of different chemical compositions is A-dependent, which 
will result in an A-dependent knee of each composition. T his feature is consistent with the 



property of the fine structures found in experimental data (lErlykin &: Wolfendald 120091 ). In 
addition, the interaction will cause a pile-up of the particles below the threshold. We expect 
this effect can contribute to the sharp knee and irregular bumps of the CR spectra. 

This paper is organized as follows. In Sec. 2 we will first go over the model describing 
CR-photon interactions. In Sec. 3 we present the calculated results and comparisons with 
observational data of the sharp knee and irregular structures of CR spectra. Finally Sec. 4 
is the conclusion. 



2. Interactions between CRs and ambient photons 

The model to explain the knee and electron/positron excesses using nuclei-photon inter- 
actions around the acceleration sources is proposed in Paper I. Here we readdress the basic 
physical picture and give some technical details. 

There are three kinds of interactions between CRs and photons: pair production, pho- 
todisintegration and photo-pion productions processes when very high energy CR nuclei 
interact with backgrou nd photons. The c ross sections for pair pr oduction and photodis- 
integration are given in iBlumenthall ( 19701 ) and IPuget et al.l (1976) respectively. The pion 



we employ an A^'^^ dependence for other nuclei with atomic number A ( IStanev et al. 



production cross section for proton is adopted from iParticle Data Group et al.l (120081). an d 



19851). 



The cross sections as functions of photon energy in the nuclei rest system fo r proton, He an d 
Fe, which are the dominant compositions for CRs around the "knee" region ( lHorandelll2003l ) . 
are shown in the left panel of Fig. [TJ Note that for proton there is no photodisintegration 
interaction. The pair production cross section is proportional to o f the nuclei, so fo r 
heavy nuclei like Fe the pair production cross section is extremely large ( IB lument hall 1 1 9 70l ) . 



To see clearly the effects of the three interaction channels, we need to know the energy 
loss of the nuclei due to each of the interactions. The energy loss rate can be written as 



dE 
dF 



kE 



kE ■ / dcos6'- 



cos 6 



de n{e) a{E, e, cos 9)c, 



where t{E) is the average interaction time of the nuclei in a radiation field with number 
density n{e), 9 is the angle between photon and nuclei momenta, c is speed of light, and 
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10" r 




Fig. 1. — Left: cross sections as functions of photon energy in the nuclei rest system for 
proton (dash-dotted), He (dashed) and Fe (sohd) respectively. For each kind of nuclei, the 
cross sections for the three interaction processes: pair production (blue), photodisintegration 
(green) and pion production (red) are shown. Right: the average relative energy loss rates 
for proton. He and Fe due to pair production, photodisintegration and pion production 
interactions in a 5000 K blackbody radiation field with respect to the nuclei energy in 
laboratory system. The labels of lines are same as in the left panel. 



K is the average fraction of energy loss in one collision, i.e. the inelasticity. We plot in 
the right panel of Fig. [T] the average relative energy loss rates with respect to the nuclei 
energy in laboratory system due to each of the interactions for proton. Helium and Iron 
respectively. The background radiation field in the calculation is specfied to be blackbody 
field with T = 5000 K. The average inelasticity k, for pair production is calculated using a 
Monte-Carlo (MC) study to average the output e+e~ pair energies for many realizations. For 
photodisintegration process the inelasticity is simply adopted as k = i/A with A the total 
number of nucleus and i the number of nucleons kicked out in one interaction. It is shown 



that i = 1 is dominant for photon energy below 30 MeV in the nuclei rest system (jPuget et al. 



1976 ). For pion production the average inelasticity is adopted as k = | ^1 -|- !!l:2_!!Li^ with 



m^r, rriA the m asses of pion meson and nuclei, and s the center-of- moment system energy 



Steckerl (119681 ). We can see from Fig. [T]that for energies lower than several PeV, the pair 
production dominates the energy losses for proton and Helium. The photodisintegration 
process of Helium become important above ~ 4 PeV. Pion production is only important for 
energies higher than several tens PeV. 

We use MC method to simulate the interactions between CR nuclei and ambient pho- 
tons. The acceleration processes and interactions are decoupled for simplicity. The CR 
nuclei are assumed to inject into the radiation field with power-law spectra and proper 
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relative abundance according to the measurements ( iHorandell l2003l : iBertaina et al.l |2008| ) . 
Note here we adopt a correction of the measured spectra to the source spectra, taking into 
account the propagation effect (see below). After interactions the CRs together with the 
e+e" products enter the interstellar environment and propagate diffusively in the Galaxy. 
For the propagation of CR nu clei we simply us e the leaky-box model with escape time 



rcsc(-R) ^ 2 X 10^ {t§v) °^ (jPutze et al.ll2009l ). The propagation of e+e is a bit com- 



plicated since the dominant effect is energy loss due to synchrotron radiation and inverse- 
Compton scattering, instead of diffusion (or to say escape). We adopt GALPROP code to 
calculate the propagation of electrons and positrons. 

The radiation field around the acceleration source is assumed to be of blackbody shape, 
with a temperature several to ten thousand Kelvin. It might be true for the young supernova 
shortly after its explosion. The intensity of the radiation field does not have to be as intense 
as the blackbody radiation. We keep the density of photons as a free parameter. The photon 
density multiplied with the interaction time gives the effective interaction probability. 

As we have stated in Paper I, the photodisintegration of heavy nuclei, especially He, 
will overproduce protons through secondary production of protons. Therefore we do not 
expect the photodisintegration to play a significant role in the interaction process. We find 
that an evolving picture with decaying and asymptotically cooling of the radiation field can 
naturally explain this requirement. Schematically speaking the photodisintegration and pion 
production have higher threshold, which requires longer time for CR nuclei to be accelerated 
to exceed the threshold. But as time goes on the temperature of the radiation field may 
become lower and lower, and the interactions of photodisintegration and pion production 
need even higher energies of CRs. Therefore the photodisintegration and pion production 
can be effectively suppressed, and only the pair production takes effect. 



3. Results 

he calculated energy spectra of individual composition, including proton. He, CNO and 
Fe, together with the observational data are shown in Fig. 2. The main parameters used in 
the MC calculation, including the relative abundances and spectral indices at the source, 
temperatures of radiation field, photon column density, and equivalent time for blackbody 
field, are compiled in Table 1. The relative abundances are adopted to reproduce the locally 
observed fluxes of each species. As same as in Paper I, only one set of parameters is enough 
to describe all spectra well. A few model parameters are slightly adjusted in order to better 
agree with the flne structures of CR spectra(see detailed explanation below). However, 
the basic features as shown in Paper I are kept unchanged. The power-law index and 
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normalization of each chemical composition are the same as that used in Paper I, except for 
CNO group we use a harde r injection spec t rum b y 0.06 and a lower flux normalization at 
1 TeV by 10% according to iBertaina et al.l ( 120081 ). which will make the flux slightly go up 
around the knee region. The temperatures of the radiation field for proton and He, which 
are the main compositions of CRs, reamin the same as in Paper I. However, we change the 
temperature for Fe from ~ 2000 K in Paper I to ~ 3500 K in this work. In Paper I we set a 
lower temperature for Fe with the purpose that the second knee around 300 — 400 PeV might 
be reproduced using the drop of Fe spectrum. If we relax this requirement we have shown 
that a wide range of the temperature for Fe can be consistent with the current data (see the 
systematic check done in Paper I). As we will see below, a higher temperature for Fe will be 
better to reproduce the fine structures of the total CR spectrum. For other nuclei which do 
not play a significant role, we use the temperature field 5500 K for 3 < Z < 25 and 3500 
K for Z > 26 respectively. It is shown that the model expectations are consistent with the 
observational data. However, since the current measurements of individual CR composition 
are actually not very good, the constraints on the model parameters are weak. 

The case for the all-particle spectra is much better than the individual composition. 
The comparison of the all-particle spectra between our theoretical prediction and the data 
is shown in Fig. 3 We can see the knee forms due to the break of He spectrum around ~ 4 
PeV. This condition would require the temperature of radiation field which He experiences is 
about 7000 K. The pile-up of He particles around the threshold point helps to better produce 
the sharp break of the knee spectrum. To investigate the break behavior more quantitatively, 
we adopt the following double power- law function to fit the CR spectrum below and above 
the knee jErlvkin fc WolfendaleHioOol ) 



I{E) = AE^ 



E 



(2) 



where 7 is the spectrum index below the knee, which changes by A7 above the knee, 5 is 



Table 1: Parameter settings in the Monte-Carlo calculation. 
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Fig. 2. — Energy spectra for proton, H e, CNQ and Fe. The sol id lines are th e MC calculated re- 



Apel et alJl2009l). AMS (I A 
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suits. Observati onal da ta are Tibet-Ill (lAmenomori et al.l 120061) . KASCAD E (lAntoni et al 



caraz et al.ll2000l : I AMS Collaboration et al.ll2002l). ATIC-2 (IPanov et al 



2001 



19931'!. CAPRICE (IBoezio et al.l l2003l). Chica go (IMueller et al 



2004), GRAPES-3 (IGupta et al.ll2008 ). JACEE (lAsakimori et al.lll998l) 



SOKOL (Ivanenko et al 



(IMueller et al.l Il991 ). EAS-TOP. H areyama (IHareyamal Il999l ). Ichi mura (llchiniura et al.l Il993l) 
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1993 ), Zatsepin (IZatsepin et al 



Swordv et al 



1995), CRN 



Juliu sson (|juliussonlll974l ). Sanriku (IKamioka et al.l 119971 ). CREAM (jYoon et al 
2OO9I ): 



200 



Ahn et al 



the sharpness parameter which describes the smoothness of the transition. The sharpness 
S (the second-order derivat ive of spectrum at the energy break point E^) is defined as 
jErlvkin fc Wolfendalel[2009[ ) 



(3) 
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Through fitting the theoretical calculated spectrum we can obtain the sharpness S = 2.0, 
while the sharpness of experimental data of Tibet air shower array is about 2.4 ± 0.8, which 
shows good consistence with the theoretical one. 




Fig. 3. — The calculated all-particle spect rum. Also shown are p roton, H e, CNO and Fe con i- 
ponents. The observational data: Tibet-II KjAmenomori et al.ll2008l ). Akeno (jNagano et al.lll984l ). 
The Normalized data are derived by corr ibing all data with a rescale based on the extrapolation of 
the direct measurements (iHorandell 120031 ) . Label "B" indicates the high energy tails of the proton 
and He spect ra above ma rker "*" can be regarded as the "component B" of the Galactic CRs as 
suggested by iHillasI (|2005l ). 



We can also note that there are high energy tails of the proton and He spectra due to 
the decrease of the pair production energy loss rate (Fig. [1]). This behavior is very similar to 
the ankle of the ultra high energy CRs generated by e"'"e~ pair production interactions with 
cosmic microwave background photons (IBerezinsky et al.ll2006l ). Coincidentall y, it may be a 
natural explanation of the "component B" of the Galactic CRs as suggested by lHillasI (120051 ) 
when studying the overall property of CR spectrum. However for the heavy nuclei, there 
are no obvious high energy tails due to large cross sections and the constant temperature of 
individual component, as shown in[T] 



Fig. m shows a more detailed comparison between theory and data, after subtracting 
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the main behavior as fitted using Eq.Q, which enables us to explore the fine structures 

of the CR spectra. The data points are c ombined using 6 experiments, Tibet Ai r Shower 

arrav JAmenomori et aliboosl). KA SKADEfjApel et al.ll2009h. ARAGATS- GAMMA(lGa rvaka et al, 



2008|), Yakut sk fllvanov et al 



2009h.M AKET-ANl Jchilingarian et allboOTh and TUNKA [56] 



according to lErlykin fc Wolfendald ( l2009l ). To avoid the uncertainty of energy calibration 
among various experiments, we fit the break energy E'^ in Eq.(|2]) for each data set and then 
normalize it to the Tibet result. It is shown in Fig. 2 that the fine structures of the theoreti- 
cal calculation is well consistent with the observational data. The position \og{E / E^) ~ 0.0 
corresponds to He, and the two bumps at \og{E / E^) ~ 0.5 and \og{E / E^) ~ 1.3 should 
correspond to CNO group and Fe respectively. The constraints on the model parameters 
are not very strong using the present observational data, e.g. the final results do not change 
significantly by varying the temperature parameters by several tens percent. 



0.15 



Tj 0.1 
o 



"m 0.05 
o 



-0.05 



Combined data from Eriykin analysis 



Monte Carlo Calculation 



knee(He) 



proton 



CNO 



— ■ — 1 ' * 



Fe 



t « 
I » 



/ 




1 




1 




1 1 


' 1 


1^ 




1 





I 

I 

I 

I 

I 

I 



-1 



-0.5 



0.5 



1.5 



LogCE/E-^)) 



Fig. 4. — The fine structures of CR spectrum: theory versus data. F means the observational or 
calculated flux of the total CR spectra, while / indicates the fitting to F according to Eg . (2). 
The data are combined results of 6 experiments according to lErlykin &: Wolfendald (|2009l ) . The 
positions \og{E/E^) ~ —0.7 and \og{E/E^) ~ 0.0 correspond to proton and He, and the two bumps 
at \og{E/E^) ~ 0.5 and \og{E/E^) ~ 1.3 correspond to CNO group and Fe respectively. 



Finally we check the spectra of e^e and position fraction of these new parameters and 
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find almost identical results as that derived in Paper I. This is reasonable because the basic 
parameters of proton and He, which produce the main part of the electrons/positrons, are 
almost unchanged. 



4. Conclusion 

In summary we use the pair production interaction model between CR nuclei and am- 
bient radiation field proposed in Paper I to explain the features of the CR spectra, including 
the sharp knee and fine structures. Results show that the spectra of CRs agree well with 
the observations. In our model, the He composition dominates around the knee at ~ 4 PeV. 
The sharp knee observed by Tibet air shower array and confirmed by more and more exper- 
iments, can be well reproduced through the pile-up of He particles. In addition, this model 
can explain the fine structures of CR spectrum through the pile-up effects of CNO group 
and Fe nuclei. As an additional result, our model can provide a natural explanation of the 
"c omponent B" problem of individual composition in tens of PeV energy range as suggested 



m 



HillasI ( 120051 ). As in Paper I, the electron/positron excesses can also be explained. 



Furthermore we would like to discuss some implications of the model parameters. We 
note that o nly one siri gle set of parameters is enough to explain the data. As discussed in 



Paper I and iHul fl2009[ ). it may indicate that the sources are "standard" which have similar 
parameters such as the temperature evolution of the radiation field, the relative abundances 
and spectral indices for individual elements, or the observed fine structures of CRs spectra are 
mainly due to one single source, either nearby or not so nearby but intensive (e.g. possibly 
the Galactic center). This work can be regarded as one progress approaching the origin of 
CRs. 
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